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Abstract 



It is demonstrated that in a two-stage scenario with elementary Poissonian emitters of parti- 
cles (colour strings) arbitrarily distributed in their number and average multiplicities, the forward- 
backward correlations are completely determined by the final distribution of forward particles. The 
observed linear form of the correlations then necessarily requires this distribution to have a negative 
binomial form. For emitters with a negative binomial distribution in the produced particles dis- 
tributed so as to give the final distribution also of the negative binomial form, the forward-backward 
correlations have an essentially non-linear form which disagrees with the experimental data. 
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1 Introduction 

Multiparticle production at high energy is currently (and successfully) described in a two-stage sce- 
nario. At the first stage a certain number of colour strings are formed, stretched between the incoming 
partons. At the second stage these strings decay into the observed secondary hadrons. In the sim- 
plest version the strings emit particles independently. However one may also introduce the interaction 
between strings in the form of their fusion and/or percolation. In any case the observed multiplicity 
distribution is a convolution of the probability distribution of the possible string configurations and 
the multiplicity distribution generated by the individual strings. In fact both are not well-known, so 
that to understand the dynamics of the multiparticle production one has to reconstruct both from 
a single distribution of the observed secondaries. To facilitate this task it is natural to use some 
other experimental information about the multiparticle production. An immediate candidate is the 
long-range correlation data, currently studied as the forward-backward correlations (FBC). They cor- 
respond to studying the average multiplicity in, say, the backward hemisphere (B) as a function of 
the event multiplicity in the forward hemisphere F. 

The most striking feature of the data is that they can be almost perfectly represented by a linear 
function [1-4] 

(B)=a + bF, (1) 

the strength of the correlation measured by the coefficient b. The dynamical reason for this simple 
behaviour is not known. In this paper we try to understand this behaviour in the two-stage scenario, 
decribed above. 

Our starting point is the assumption that individual emitters (colour strings) have a Poisson 
distribution in the number of produced particles. This assumption is based on the idea that a string 
is an homogeneous extended object in the rapidity space. Its fundamental property is that parts of 
the string occupying different regions in the rapidity space do not interfere (except at small rapidity 
distances). Then the particle distributions coming from parts of the string should have basically the 
same form which should also be similar to the overall particle distribution (except for the average 
multiplicities, which shoud sum together into a total one). In terms of the corresponding generating 
functions (GF) this means that the functional dependence of the product of GF corresponding to 
different parts of the string should be the same and coincide with the functional dependence of the 
total GF for the string. As we shall see, this condition is fulfilled by the exponential dependence of 
the Poisson GF (see Eq. (22)). It is difficult to imagine any other probability distribution with the 
same properties. 

Our main result is that with the individual emitters having the Poisson multiplicity distribution, 
the FBC are completely determined by the final multiplicity distribution in the forward rapidity 
window P(F) and do not depend explicitly on the distribution of strings themselves (and consequently 
on their dynamics). If, in accordance with the experimental data, one further requires that (B) be 
a linear function of F, then one finds that the corresponding P(F) must be a negative binomial 
distribution (NBD), also in agreement with the experimental data. Thus our explanation of the linear 
dependence (1) is that it comes from a Poissonian form of the elemental distribution, which the 
distribution in the elemental emitters transforms into a NBD. As we shall see such transformation 
is only possible with different kinds of strings, including those with very small multiplicities. Such 
strings are indeed present in current scenarios, either as "short" strings, extended over comparatively 
small rapidity intervals, in the independent string picture, or strings with small transverse area, which 
appear in the percolation scenario. 

This result is not at all trivial, as one can conclude from an alternative picture in which both the 
strings and the final distribution are negative binomial. In this case we find that (B) is not a linear 
function on F, but rather a function which grows with F with a falling derivative, so that at high F 
it either saturates or grows only logarithmically. 

To compare our predictions with (scarce) experimental data we take into account that most (if 
not all) particles are created as resonances. So in fact our numbers F and B refer not to the observed 
particles but only to the directly formed ones. As in most other approaches we then have to introduce 
an average number r of the observed particles per single directly produced one (resonance, or cluster). 
The number of observed particles in the forward or backward rapidity windows will then be iif = rF 
and iib = rB respectively. Of course, from (1) also a linear relation will follow between (ns) and np 
with the same parameter b. Using the experimental data on the forward multiplicity distribution in 
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[5] and comparing them with the data on the parameter b [4] we find that r does not seem to depend 
on energy and its value is around 3. 



2 Formalism 

Let us consider a set of N identical elemental emitters distributed with a probability w(N). We denote 
p(F, B) the probability to find F(B) particles in the forward (backward) hemispheres from a single 
emitter. The overall distribution in F, B coming from any number of emitters is given by 

JV 

P(F,B) =^>(A0 £ S F ^ Fz S B ^ Bz l[p(F l ,B l ). (2) 

N Fi,B t i=l 

The overall averages are 

(A)=^A(F,B)P(F,B). (3) 

F,B 

The conditional probability to find B particles in the backward hemisphere with their number F in 
the forward hemisphere fixed is 

rm - (4) 

Here P{F) is the (unconditional) probability to observe F particles in the forward hemisphere 

P{F) = Y J P{F,B). (5) 

B 

According to (4) and (5) the average with a given F is given by the ratio 

_ j: B A(F,B)P(F,B) 

{A)f - EbWb) • (6) 

It is convenient to pass to the formalism of generating functions (GF). One represents the Kronecker 
symbols entering (2) as 



where the integration contour in the complex x plane goes around the origin. Using this formula one 
writes for the probability (5) 

™°>-J£;*-'J£?-"Z •"<"» S ( E »**«.*>) • e» 

Each sum inside the product over i = 1, ...N gives a GF for the distribution p(F, B) from a single 
emitter: 

g(x,y) = Y,x F y B p(F,B). (9) 

F,B 

So we get 

where we introduced the overall GF 

G(x,y) = '£ i w(N)g N (x,y). (11) 

N 

The averages with a fixed F can be readily expressed in terms of G. We find 
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To do the sum over B we use 



?4W--A(*-^)iy*-*(*-.£)^r 



(13) 



The integration by parts in y then gives 



-Lmmwrn- j J ' 53*^ (*»|;) <««.»>■ <14 » 



B 

Taking finally the residue at y = 1 we find 



s 

In partcular we get 



and 



^(F,B)P(F,B)=I ^- F [A(F,y^)G( X ,y)]^ (15) 

P(F) = £ P(F, B) = J ^-x~ F G{x, 1) (16) 
^BP(F,B)= r^L X -*G' y ( X ,l). (17) 



This allows to write the average number of particles in the backward hemisphere, given their number 
in the forward hemisphere, as 

/ dxx- F - x G' y {x,\) 
{B)F = Jdxx-^G(x,l)- (18) 

This formalism can be generalized in a straightforward manner to the case when the emitters are not 
identical. Call a configuration C a particular number Nq of emitters enumerated by i — 1,2, ....N c 
all of them (or part of them) different. Let the probability to find F(B) particles in the forward 
(backward) hemisphere from the z-th emitter be Pi(F, B). Define corresponding GF gi(x, y) as in (9). 
Then for the overall GF one finds an evident generalization of (11) 

N c 

G(x,y) = J2w(C)]Jg i (x,y), (19) 

C i=l 

where w(C) is the probability to find a given configuration C and summation goes over all possible 
configurations. 

3 Poisson emitters 

In the following we assume that individual emitters have no FBC, which implies that both the prob- 
ability and GF for them factorize 

p(F,B)=p(F)p(B), g(x,y) = g(x)g(y). (20) 

As to the concrete form of the distributions p{F) and p(B), we assume them identical and Poissonian: 

p(F) = e-^. (21) 

The corresponding GF is an exponential: 

g(x,y) = g(x)g(y), g(x) = e*^. (22) 

The reason to assume the Poissonian distribution for a single emitter has been explained in the 
Introduction: it naturally goes with the notion of the string as a homogeneous extended object in the 
rapidity space. 
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Passing to the overall GF G, we see that it does not generally factorize. With (22) we find 

G(x, y) = Y, w{N)e N < x+ y- 2) . (23) 

JV 

The only case when it does factorize corresponds to a fixed number of emitters w(N) = 5n,n - 

The special property of the Poisson distribution for the elemental emitters shows itself in that the 
overall GF results depending only on the sum x + y. As a consequence we find 

G' y (x,l) = G' x (x,l). (24) 
Putting this into (17) and integrating by parts in x we transform the numerator in (18) as 

(F + 1) / ^-x~ F - x G{x, 1) = (F + 1)P(F + 1). (25) 

As a result (B) p becomes expressed entirely in terms of the distribution in the number of the forward 
particles: 

_ (F+1)P(F + 1) 
{B)F P{F) ■ ( } 

Before we go further, note that this result holds not only for identical emitters, but also for different 
but Poissonian emitters. In fact the key point is that for an individual Poissonian emitter 

Pi (x,y)=e a ^ + y-V = Pi (x + y) (27) 

depends only on the sum x + y. Summation over configurations in (19) will not change this result, so 
that the crucial identity (24) and the following equations will be valid. 

Eq. (26) allows to find (B) once the particle distribution P(F) in the forward rapidity window is 
known. A more detailed form of Eq. (26) can be obtained in the approximation of the KNO scaling, 
when (B) can be expressed through the scaling function (see Appendix 1.). 

4 Linear form of (B)p 

As we have shown, under the assumption that elemental emitters are Poissonian (B) p can be expressed 
entirely in terms of the distribution in the number of the forward particles by Eq. (26). Let us assume 
that the dependence of (B) p is linear and given by (1). Then from (26) we obtain a recurrent relation 

P(F + 1) = ^^-P(F) = (b+ j^j P(F), (28) 

where c = a — b. As we demonstrate below, it allows to completely determine P(F) (see also [6]). 
We put 

P(F) = b F Q(F). (29) 



Then we find a recurrency for Q(F) 



F + l 

where X = c/b = a/b — 1. Recurrency (|30[) is trivially solved by 



Q(F + l)=ll + -±-)Q(F). 



lnQ(F) =lnQ(0) + ^ln(l + -j . (31) 



The sum over n can be expressed via the Gamma- function. Indeed we find that if 

F 



n— 1 ^ ' 
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then its derivative in A 

I i 

S'(A) = ^— -=V0F + A + 1)-^(A + 1). (33) 



, n + A 



Integrating this over A in the interval [0, A] we get 



r(f + A + i) 

* = ln r ( F + i)r(A + i) (34) 



This finally gives for the distribution P(F): 



P(F) = P(0)e Fhlfc r(F+ ,M 1} r , (35) 

y ' w r(F + i)r(A + i) v ' 

where P(0) is to be determined from the normalization condition. Evidently one should have b < 1 
for normalizability. 

One observes that ( |35| ) is nothing but a negative binomial distribution with the average value (F) 
and parameter k determined by 

b = TFjTk> * = * + ! = £ ( 3fi ) 

Thus the NBD is singled out from all other distributions by the fact that if the elemental emitters are 
Poissonian it leads to a linear dependence of (B) on F. However to give this result a physical meaning 
we still have to demonstrate that it is possible to produce a NBD distribution by a convolution of 
Poisson distributions with a physically admissible distribution of emitters. This problem is discussed 
in the next section. 



5 Negative binomial distribution from Poisson distributions 

We start by showing that if all the emitters are identical and Poissonian, one cannot produce a NBD 
distribution of final particles with any physically admissible distribution of emitters. Indeed let P{F) 
have a NBD form. The corresponding GF is then 

G N bd(z) = {a-bz)- k , a- 6 = 1. (37) 

The distribution itself can be extracted by developing this function in a power series: 



Pnbd(F) = a~ C£_ 1+F 1 — 1 , (38) 

where we define the binomial coefficient as 

F T(F + fc) 

^k-l + F- T{k)T{F + 1} W 

The average value of the multiplicity is given by 

(F) = k(a - 1), (40) 

so that the NBD is fully determined by two parameters: the average multiplicity and k. 

Now it is trivial to show that if the GF are known for both the individual emitter and for the overall 
distribution, the GF for the distribution of emitters is determined uniquely (provided all emitters are 
identical). Indeed, considering again the distribution of only forward particles, we have seen that the 
overall GF is given by 

G(z)=^(N)g N (z), (41) 

N 
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where g(z) is the GF for the individual emitter. Introducing the GF for the distribution of emitters 

by 

H(z) = J2 w ( N ) zN > ( 42 ) 

N 

we observe that 

G(z)=H(g(z)). (43) 

Assume that we know both G(z) and g(z). We solve g — g(z) for z to determine z = z(g). Putting 
this into G{z) we find the unknown function H{z): 

H(g) = G(z(g)). (44) 

As we see, the GF for the emitter distribution is determined in a unique manner. However it may 
be or may not be physical, depending on its behaviour at small values of argument. For a physically 
valid H (z) it should be analytic at small z with all the derivatives positive. Otherwise it cannot be 
interpreted as a GF for a physically sensible distribution. 

Let us see what this means for the case when the final distribution is negative binomial and the 
individual emitter is Poissonian. Then the GF for an individual emitter is given by (22) Inverting we 
find 

2 = 1 + — ln^. (45) 
a 

So with the final GF given by (37), the GF for the emitter distribution turns out to be 

r / 1 \ i —k r 1 i —k 

H(g)= o-6(l + - lno) = 1-6-lno . (46) 
L \ a J J L a J 

This function is not analytic at g = 0, all its derivatives being infinite at the origin. So it cannot be 
interpreted as a GF for some physically sensible distribution of emitters. 

Thus there does not exist a distribution of identical Poissonian emitters which gives an overall NBD. 
However this is not true if the emitters are allowed to be different and, in particular, continuously 
distributed in their average values cti. Indeed, for different Poissonian emitters, we have for the GF 
for the forward hemisphere 

N C 

G(z) = W ( C ) J] e"'* 2 - 1 ' = ™(c)e a(C)( " 1) , (47) 
c i=\ c 

where 

N c 

a(C) = J^a i (48) 

and has a meaning of the average forward multiplicity for a given configuration C . The quantity to 
be averaged over configurations in (47) depends only on a(C). Assuming it to vary continuously, one 
can rewrite (47) as 

/>oo 

G(z)= / daw{a)e a{z ~ 1 \ (49) 
Jo 

where w(a) > now has the meaning of the probability distribution for the average forward multi- 
plicity of the emitters, normalized accordingly 

poo 

/ daw(a) = 1. (50) 
Jo 

Now it is trivial to find a form of w(a) which leads to the NBD. In fact let us take 

W) 



k 

w (a) = ^—a k - 1 e ra . (51) 



Then (47) gives 



GW=( 7TT — ) ={a-bz)-\ (52) 
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where 

a=l + l/r, b=l/r. 

So we have indeed obtaind the NBD (37) . From the derivation it can be seen that the crucial point in 
obtaining a GF which behaves as a meromorphic function, instead of the entire Poissonian, has been 
the integration over small valuses of a down to a = 0. Physically this implies existence of emitters 
(strings) with very small average multiplicities. 

Note that a construction similar to (49) with (51) was used in the VENUS model [7] to obtain a 
realistic NBD for p-A spectra. 

To conclude, in the picture where the Poissonian emitters arc different and may have very small 
average multiplicities, it is possible to choose their distribution to have the overall NBD. In this 
picture the FBC have an exact linear dependence as in the empirical formula (1) with the parameter 
b determined by (36) 

6 Negative binomial distribution from negative binomial dis- 
tributions 

Although the results of the preceding section show that one can assume the elemental emitters to be 
Poissonian and end up with the NBD for the finally produced particles, provided the emitters are 
different and extend to very small average multiplicities, it presents some interest to see what happens 
in a more conservative picture, in which all elemental emitters are identical. As we have seen, to 
have the final NBD the emitters cannot be Poissonian in this case. So the first problem is to find an 
adequate distribution for a single emitter which would lead to the NBD with an approprate physically 
meaninful emitter distribution. This problem can evidently have many solutions. One of the simplest 
was found by C.Iso and K.Mori [8]. In the following we use their results to study the long-range 
correlations which follow from their form of the distributions. We expect that the conclusion will not 
depend too heavily on this particular choice. In any case they will serve as a simple example to be 
compared with the Poissonian emitters. The solution of C.Iso and K.Mori is to take the GF for the 
individual emitter essentially also of the NBD form 

9 {z) = - d -\{a-bz)-v, (53) 

with the requirements: 

y>0, d = e-l<0, e = a~ y . (54) 

The latter condition follows from the requirement that the emitter should emit at least one particle, 
i.e. g(0) = 0. The average multiplicity from a single emitter corresponding to (53) is given by (F)/xd 
, where (F) is the overall average corresponding to the NBD. With (53) one easily shows that the GF 
for the emitter distribution has a form 

H(z) = (e-dz)- x , (55) 

where x < is determined by the final BND parameter k: k = — xy. It turns out that (55) corresponds 
to a physically meaningful distribution so long as — x > 2. The solution of C.Iso and K.Mori contains 
only one more parameter y (or x) in addition to the two parameters of the final NBD. 

Passing to the FBC, with a non-Poissonian emitters our simple formula (26) does not hold and we 
have to explicitly calculate them from (18). We have to calculate 

G' u (z 1 l)=\Y J ™{N)g N (z)Ng N ~ 1 (u)h'(u)\ . (56) 

j u,=i 

Using the explicit form of g{u), Eq. (57) we find 

9'{l)=\-^{a-bu)-y-A = * (57) 
L a J u=i a 



We are left with 



G'u(z, 1) = E ™(N)Ng N (z) = -%g± £ w(N)g N = - ^(e - dz)~*. (58) 
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In the last equation we used (55). Taking the derivatives and putting (37) for g(z) we finally obtain 



hk r 

G' u {z, 1) = [e(a - bz)- k+y - (a - bz)- k 



(59) 



Both terms are the GF for the NBD, with parameters k — y and k respectively. So projecting their 
F-th term in the development in the power series in z will give 

eP N BD(F, k-y)- P NBD (F, k). 

According to (18) to obtain (B) we have to divide this by Pnbd {F, k) . This brings us to our final 
result 



(B) 



bk 
1 



P N BD(F,k-y) 



Pnbd{F, k) 

or using relations between e, d and the overall average (F) 



- 1 



(B) 



(F) 



1 - a~y 



1 



- v Pnbd(F, k-y) 
Pnbd (F, k) 



(60) 



(61) 



where 



k 



+ 1 



and y is the only extra parameter introduced in the solution. Using the explicit expression for the 
NBD distribution we find 



w \-a-y 



1 



T{k) T(F + k-y) 



T(k-y) T(F + k) 



(62) 



At large values of F the negative second term is vanishing as (F + k)~ v . So for y > (B) grows 
with F tending to a constant value, which is somewhat greater than its overall average. In the limiting 
case y = one gets 



(B) 



(F) 



ip{F + k) -ip(k) 



(63) 



So (B) continues to grow indefinitely with F, but only logarithmicallly. Thus, comparing with the 
Poisson emitters case, we find the same tendency in the behaviour of (B): it grows with F. However 
now the growth is either stopped at high F or severely moderated. 



7 Comparison to the experiment. Discussion. 

Passing to the comparison with the latest experimental data [4], we first have to stress that the 
observed values of b strongly depend on the rapidity windows chosen. They diminish substantially 
with the gap between the forward and backward windows. This is thought to be an effect of the 
remaining short range correlations. However the fact that b continues to considerably go down when 
the gap is raised from 4.37 to 5.43 makes this explanation dubious. In any case there is little reason 
to think that the FBC are universal for different choices of the rapidity windows, so that the analysis 
has to be made for each particular choice separately. To get rid of the short range correlations the gap 
should be taken large enough. From this point of view the data which involve the full forward and 
backward spaces are evidently not suitable. On the other hand our results are presented in terms of 
the particle distribution in the forward window. Such data are rather scarce in general and still more 
scarce if we want to have them exactly for the windows which have been used for the determination 
of the FBC. Thus, comparing the data from [4] and [5] we find only one pseudorapidity window 

2.18 < \r]\ < 3.25 (64) 

in which the parameter b has been measured [4] and which is covered by two pseudorapidity windows 
of nearly the same length 

2.5 < \r)\ < 3.5, and 2.0 < \t}\ < 3.0 (65) 

where the particle distribution has been studied [5] . The data on the FBC in the window agree with 
the form (1) with the parameter b given in the Table 1 for various cm, energies. As to the particle 
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distributions in the two rapidity windows (65), it is claimed to be well described bu the NBD with 
parameters listed in Table 2 for two cm. energies 200 and 900 GeV. As we see neither the windows 
nor the energies do not match in the two sets of data. In the following we use the values of the NBD 
parameters for the interval (64) as the average for the two intervals (65) and also use the parameters 
in Table 2 for energies 200 and 900 to desribe b at 300 and 1000 GeV respectively. 

Turning to our results we first stress that the form of (B) which follows from the NBD for individual 
strings (62) or (63) is in striking contradiction with the data. In Figure we illustrate this behaviour for 
two possible values of the parameter y = and 0.5 and compare it with the straight line corresponding 
to b = 0.252 and a = 1.2 ([4], with an appropriate reduction of a for the interval (64)). The number r 
of observed particles per directly produced was taken unity. However changing it does not improve the 
situation at all. So our first conclusion is that the scenario of strings with the NBD, which produce 
the overall NBD, does not work. 

In our basic picture with Poissonian strings the linear form (1) holds exactly. The value of the 
parameter b is given by Eq. (36) where the average (F) refers to the directly produced particles. 
Passing to the observed particles we find 

{HF) (66) 



(n F 



rk ' 



where r is the number of observed particle per one directly produced. The data shown in Tables 1 
and 2 allow to extract the value of r for the rapidity window (64) and two cm. energies 200 and 900 
GeV. We find 

r = 3.22 ±0.68, y/s — 200 GeV, 

r = 2.63 ± 0.52, s/1 = 900 GeV. (67) 

Note the large errors, which come from the errors in the experimental b and fitted k. From these 
values it is difficult to extract any definite information as to the behaviour of r with energy. They are 
compatible with r independent of energy altogether. 

If we forget about the necessity to cut the central rapidity window and consider the whole phase 
space both in the forward and backward hemispheres, we get values of r with less errors. They are 
presented in the first column of Table 3 for various energies. They also seem to be independent of 
energy. 

Turning to the literature, we find that the FBC have been mostly discussed for the whole phase 
space in terms of the assumption by Chou and Yang [9] about the form of the multiplicity distributions 
in the total number of particles n and in the difference z — np — ub- On purely phenomenologal 
grounds Chou and Yang assumed that this distribution is a product of the NBD in n and a binomial 
in z, the latter proportional to C™jJ r , where r is the number of particles per cluster. From this 
assumption, it was found in [10] by numerical calculation that the resulting behaviour of (B) as a 
function of F was close to linear. In this approach the parameter b can be expressed through the total 
multiplicity as 

b= { ;\ +k - r . (68) 

(n) + k + r 

Comparing this formula with the experimentally observed NBD parameters the values of r were found 
for various energies in [4]. They are presented in the second column of Table 3. As we see they are 
smaller than ours and seem to grow with energy. 

Commenting on these results we have first to mention that our scenario with Poissonian strings 
automatically leads to the distribution in the difference z — F — B of the binomial type (see Appendix 
2.). So in this way the formal assumption of Chou- Yang receives a concrete realization in our scheme. 
Since the linear dependence (1) is exact in our approach, we have little doubt that it should also hold 
in [10], which the numerical calculations made there confirm. However there remain a quantitative 
difference in the two approaches. One observes that for r different from unity the final formulas for 
b (36) and (68) remain different if one considers the whole phase space and takes into account that 
(n) = 2(F). Only for r = 1 they coincide. As a result, our values for r are systematically larger 
than those obtained in [4] from the essentially Chou- Yang approach (cf. the two columns in the Table 
3). On the other hands our values do not essentially change with energy, whereas those of [4] grow, 
so that the difference between them seems to be disappearing at larger energies. The exact origin 
of this difference is easy to trace. It is related to the way in which the distribution in clusters is 
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transformed into the one for the observed particles. In both approaches this transformation is made 
in an approximate manner. In the purely phenomcnological approach of Chou-Yang the distribution 
in the total number of observed particles is taken to be the NBD just on experimental grounds. 
However it does not evidently correspond to the idea that particles are formed from clusters. In 
fact if we assume that all particles are formed from clusters which decay in exactly r particles, their 
distribution cannot be an NBD in n, but rather an NBD in n/r. In this case one will only find a 
non-zero probability to see r, 2r, 3r... observed particles. In our approach the distribution in n is taken 
to be an NBD in n/r suitably continued to non- integer values of the argument. This is evidently also 
only an approximation. It is trivial to find the exact distribution in n which corresponds to cluster 
formation and decay in our approach (see Appendix 2.). Understandably it is not an NBD in the 
general case, passing into it only when the number of particles into which a cluster may go is fixed. 

In conclusion wc think that due to the introduction of a certain dynamical mechanism for particle 
production via interacting strings our approach allows not only to understand the reason for the formal 
assumption made in [9] and exploited in [4,10], but also presents some possibilities to improve the 
description of the FBC. In particular concrete models for string fusion and percolation introduce well 
defined distributions of strings and so allow, in principle, to predict the final form of the observed 
particle distribution. Unfortunately quantitative realization of this idea necessarily requires taking 
into account energy-momentum conservation at each step, which can only be achieved numerically, in 
a well developed Monte Carlo algorithm. 
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9 Appendix 1. KNO scaling functions 

Using the KNO scaling for P(F) we present it as 



(69) 



with some smooth scaling function ip(x). Then 



P(F + 1) = (F)- V = P(F) + (F)-y (^y) • (70) 

Putting this into (26) we obtain 

(B)F=F + l + ^±±[ln4,(x)}' x=F/{F) . (71) 

This formula is generally quite different from (1). Its exact form depends on the form of the KNO 
scaling function. Let us assume 

i/>(x) = Aer^- 1 ^. (72) 

Then 

\}nip(x)]' = -sign(x - l)a/3\x - l| a_1 , (73) 

so that we finally get 

(B) F = F + 1 + sign(F - (F»a/?(| JL p - | ^ |") . (74) 

As one observes, as a rule, the F dependence is much more compliciated than the linear one assumed 
in (1). 



12 



10 Appendix 2. The distribution in z and production of clus- 
ters 

We first show that the Poisson distribution of particles produced from the individual emitter automat- 
ically leads to a binomial distribution in z at fixed total number of produced particles. (It is trivial to 
show that also the inverse is true: the binomial distribution in z means that the elemental emitters are 
Poissonian). We again use the formalism of GF. In this language transition from numbers of particles 
F and B to n = F + B and z = F — b corresponds to the transition from variables x and y to u and 
v determined by 

x = uv, y — u/v. (75) 
The GF G(u, v) is related to the distribution P(n, z) via 

n 

%«)=EE«V. (76) 

n— z=—n 

Now assume that the elementary emitters are Poissonian so that G(x, y) = G(x + y)- We present it as 



G(x + y)=J2 < n )( x + y) N = J2 a(n)(«/2)> + W - £ a(n)(«/2) n £ 



2m 
m=0 



= 5>(n)(u/2)» £ C^V. (77) 

n z=—n 

From this expansion we conclude that 

P(n, z) = (1/2)" 'a(n)C^' 2 = (l/2)»a(n)C^ (78) 
and the overall distribution in n is 

n 

P(n)= ]T P(n,z) = a(n). (79) 

z= — n 

Dividing (78) by (79) we get the distribution in z for fixed n , which results proportional to C% ■ 

Now we discuss the correct treatment of resonance formation. Let the elemental emitter (string) 
produce two types of resonances 1 and 2 with a probability p(ni,n 2 ). Let the resonances decay into 
r\ and r 2 observed particles respectively. The final distribution of the observed particles will be given 
by an evident generalization of Eq. (2): 

P(n) = 5>(tf) £ 6 n _ Ei{rinii+r2n2i) .l[p(nu,n 2i ) (80) 
The GF for the distribution in n will be given by 

If we introduce the GF for a single emitter 

ff (tti,u 2 )=5]p(ni ) n2K 1 ^ ) (82) 

then Eq. (81) can be rewritten as 

G(u)=J2^(N)g N (u r \u^). (83) 

Now assume that the two resonances are produced independently and the distribution in both is 
Poissonian: 

g(u 1 ,u 2 ) = e «i("i-i)+«2(«2-i) ! ( 84 ) 



G(u) = P{n)u n = £ w(N) ( p{n 1 ,n 2 )u rini+r2n A . (81) 
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where 0^2 give the average number of each sort of resonance produced by a single emitter. Also 
assume that the emitters are distributed like in (49) with the ratio o^/ai fixed. Then in complete 
analogy with (52) we shall get an overall GF in the form 

G(u) = (a- b lU ri - 6 2 w r2 ) k (85) 

with a — b\ — 62 = 1- Evidently this is not a NBD. It passes into a NBD only in two limiting cases: 
b 2 = when it gives a NBD in the observed particles and 61 = when it is a NBD in resonances and 
n takes on values which are multiples of r. 
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12 Figure caption 

(B) as a function of F with the NBD for individual emitters. The lower curve corresponds to y = 0.5 
in Eq. (62), the upper one corresponds to y = (Eq. (63)). The straight line shows the experimental 
behaviour at cm. energy 200 GeV. 

13 Tables 

Table 1. Experimental data on b for the rapidity window (64) [4] 



Energy (GeV) 


300 


546 


1000 


1800 


b 


0.252 ±0.044 


0.376 ±0.033 


0.381 ±0.0021 


0.501 ±0.008 



Table 2. The NBD parameters for the rapdity windows (65) [5] 



Energy (GeV) 


^-interval 


(F) 


k 


200 


2.0-3.0 


2.40 ±0.03 


2.01 ±0.07 ±0.10 


200 


2.5-3.5 


2.14 ±0.03 


2.2 ±0.1 ± 0.1 


900 


2.0-3.0 


3.72 ±0.05 


2.2 ±0.1 ±0.1 


900 


2.5-3.5 


3.41 ±0.05 


2.2 ±0.1 ±0.1 



Table 3. Number r of particles per produced cluster 



Energy (GeV) 


r, Eq. (66) 


r, Eq (68) [4] 


300 


2.86 ±0.29 


2.27 ±0.24 


546 


2.73 ±0.26 


2.31 ±0.19 


1000 


3.11 ±0.30 


2.65 ±0.25 


1800 


2.96 ±0.21 


2.64 ±0.16 



